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1. Introduction and Summary 

Energy dissipation of a heavy quark moving through a hot plasma is both theoretically inter- 
esting and experimentally relevant [1-17]. A high energy particle moving through a thermal 
medium is an example of a non-equilibrium dissipative system. The particle will lose en- 
ergy to the surrounding medium, leading to an effective viscous drag on the motion of the 
particle. In a weakly coupled quark-gluon plasma, the dominant energy loss mechanisms 
are two-body collisions with thermal quarks and gluons, and gluon bremsstrahlung (see, for 
example, Ref. [11]). Which mechanism dominates depends on the rapidity of the quark. 
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1.1 Heavy ion collisions 

Most previous work on the rate of energy loss by a charged particle moving in a plasma is based 
on perturbative weak-coupling approximations [1-12,14], but one would like to understand 
the dynamics in strongly coupled plasmas. The specific question of the energy loss rate 
of a moving quark in a strongly coupled non-Abelian plasma is of more than theoretical 
interest. At RHIC, collisions of gold nuclei at 200 GeV per nucleon are believed to produce 
a quark-gluon plasma which, throughout most of the collision, should be viewed as strongly 
coupled [18,19]. For the early portion of the collision (but after apparent thermalization) a 
temperature T « 250 MeV is inferred, with a strong coupling a s on this scale of perhaps 0.5. 

Charm quarks (as observed through the production of D mesons) provide several im- 
portant observables. One involves the elliptic flow, denoted V2{pt), which is a measure of 
the azimuthal anisotropy of produced hadrons with respect to the reaction plane. The mea- 
surement of a large elliptic flow for light hadrons is one of the significant pieces of evidence 
supporting the claim that the quark-gluon plasma produced in RHIC collisions responds like 
a nearly ideal fluid with a small mean free path [20-28]. Since the mass of a charm quark, 
m w 1.4 GeV, is large compared to the temperature, one naively expects charm quarks to 
equilibrate more slowly than light quarks. Because elliptic flow is primarily generated early 
in the collision, slow thermalization of charm quarks should imply diminished elliptic flow for 
charmed hadrons. The extent of the delay in thermalization, and the resulting suppression 
of elliptic flow, depends on the charm quark energy loss rate [11]. 

Another detectable effect sensitive to the rate of energy loss of quarks moving through 
the plasma is jet quenching. Within the ball of plasma formed by the collision of two large 
nuclei, a quark (or antiquark) from a qq pair produced near the center of the ball is less likely 
to reach the edge with sufficient energy to form a detectable jet (after hadronization) than a 
quark from a qq pair formed near the surface of the expanding plasma. But if one quark from 
a pair created near the surface escapes and forms a jet, then the other quark, recoiling in 
the opposite direction, will typically have to travel much farther through the plasma before 
it can escape. If the rate of energy loss to the plasma is sufficiently large, then one should, 
and in fact does, see a suppression of back-to-back jets (relative to pp collisions). A related 
quantity also sensitive to this effect is the suppression factor J?aa (pt)j which is the ratio of 
the D meson spectrum in Au-Au collisions to that in pp collisions. 

1.2 AA = 4 supersymmetric Yang-Mills theory 

In this paper we present a calculation of the energy loss rate for quarks moving through a 
plasma of Af = 4 supersymmetric SU(N C ) Yang-Mills theory (SYM), in the limit of large 't 
Hooft coupling, A = gy M N c S> 1, and a large number of colors, N c — ► oo. The quarks whose 
dynamics we will study are fundamental representation particles introduced into M = 4 SYM 
by adding an M = 2 hypermultiplet with arbitrary mass. 1 In the large N c limit, fundamental 

1 More explicitly, this means adding a Dirac fermion and 2 complex scalars, all in the fundamental repre- 
sentation, with a common mass and Yukawa interactions which preserve N = 2 supersymmetry. 
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representation fields have negligible influence on bulk properties of the plasma. One may view 
the quarks as test particles which serve as probes of dynamical processes in the background 
N = 4 plasma. 

The reason for studying M = 4 super- Yang- Mills is simple — it is easier than QCD. There 
are no good approximation techniques which are generally applicable to real-time dynamical 
processes in strongly coupled quantum field theories. Thermal relaxation or equilibration 
rates, such as the energy loss rate of a moving heavy quark, cannot be extracted directly 
from Euclidean correlation functions and hence are not accessible in Monte Carlo lattice sim- 
ulations. 2 But for the specific case of Af = 4 SU(N C ) supersymmetric Yang-Mills theory, the 
AdS/CFT conjecture (or gauge/string duality) states that this theory is exactly equivalent 
to type IIB string theory in an AdS§ x S 5 gravitational background, where AdS$ is five 
dimensional anti-de Sitter space and S 5 is a five dimensional sphere [30-32]. 3 At large iV c 
and large A, the string theory can be approximated by classical type IIB supergravity. This 
approximation allows completely nonperturbative calculations in the quantum field theory to 
be mapped into problems in classical general relativity. In this context, raising the tempera- 
ture of the gauge theory corresponds to introducing a black hole (or more precisely, a black 
brane) into the center of AdS§ [33]. According to the AdS/CFT dictionary, the Hawking 
temperature of the black hole becomes the temperature of the gauge theory. 

At zero temperature, the properties of J\f = 4 supersymmetric Yang-Mills theory are com- 
pletely different from QCD. AA = 4 SYM is a conformal theory with no particle spectrum or 
S-matrix, while QCD is a confining theory with a sensible particle interpretation. But at 
non-zero temperatures (and sufficiently high temperatures in the case of QCD), both theories 
describe hot, non-Abelian plasmas with Debye screening, finite spatial correlation lengths, and 
qualitatively similar hydrodynamic behavior [34]. The major difference is that all excitations 
in AA = 4 SYM plasma (gluons, fermions, and scalars) are in the adjoint representation, while 
hot QCD plasma only has adjoint gluons and fundamental representation quarks. There are 
a variety of reasons to think that many properties of strongly coupled non-Abelian plasmas 
may be insensitive to details of the plasma composition or the precise interaction strength. 
In N = 4 SYM, bulk thermodynamic quantities such as the pressure, energy or entropy den- 
sities, as well as transport coefficients such as shear viscosity, have finite limits as the 't Hooft 
coupling A — > oo. The pressure divided by the free Stefan-Boltzmann limit (which effectively 
just counts the number of degrees of freedom) in N = 4 SYM is remarkably close to the cor- 
responding ratio in QCD at temperatures of a few times T c where it is strongly coupled [35]. 
The dimensionless ratio of viscosity divided by entropy density equals 1/4-7T in AA = 4 SYM, 
as well as in all other theories with gravity duals [36, 37] in the strong 't Hooft coupling 
limit. And this value, which is lower than any weakly coupled theory or known material 
substance [34], is in good agreement with hydrodynamic modeling of RHIC collisions [18,19]. 
These features have led some authors to speculate about "universal" properties of strongly 

2 See, however, Ref. [29] for a recent effort to extract an estimate of the damping rate by fitting a 
parametrized model of the spectral density to lattice data for the Euclidean current-current correlator. 
3 This conjecture is unproven, but is supported by a very large body of evidence. We assume its validity. 
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coupled plasmas. 

In the dual gravitational description, a fundamental representation hypermultiplet cor- 
responds to the addition of a D7-brane to the black hole geometry. This D7-brane wraps an 
S 3 C S 5 and wraps all of the Schwarzschild-AdS geometry down to a minimal radius (which 
is dual to the mass of the quark) [38]. The addition breaks the amount of supersymmetry 
in the theory down to N = 2. According to the AdS/CFT dictionary, an open string whose 
endpoints lie on the D7-brane is a meson, with the endpoints of the string representing the 
quarks. At non-zero temperature, one can also have open strings which stretch from the 
D7-brane down to the black hole horizon. The existence of such solutions reflects the fact 
that N = 4 super- Yang- Mills theory, at any non-zero temperature, is a deconfined plasma in 
which test quarks and antiquarks are not bound by a confining potential. We will extract 
the energy loss rate of moving quarks in this gauge theory by studying the behavior of the 
endpoints of both types of such open string configurations. To our knowledge, this is the first 
quantitative, nonperturbative calculation of the energy loss rate of a moving massive quark 
in any strongly coupled quantum field theory. 4 

1.3 A toy model, and the plan of attack 

The following toy model helps to clarify some of the issues which will arise in our analysis of 
heavy quark damping. Consider a particle with momentum p moving in a viscous medium 
and subject to a driving force /. Its equation of motion may be modeled as 

P = -VP + f, (i- 1 ) 

where ji is the damping rate (or friction coefficient). To infer information about fi from 
motion of the particle, it is useful to consider two different situations. First, if one examines 
steady state behavior under a constant driving force, then p = implies that p = f / fi. If the 
particle has an (effective) mass m and its motion is non-relativistic, so that p = mv, then 
the limiting speed v = //(m/x). Hence, a measurement of the steady state speed for a known 
driving force determines the combination m/x, but not fi alone. 

Second, if the driving force / = 0, then a non-zero initial momentum will relax expo- 
nentially with a decay rate of pit) = p(0) e _/lt . If momentum is proportional to velocity, 
then the speed of the particle will show the same exponential relaxation. A measurement of 
p/p, or v/v, will thus determine the damping rate [i. The important point is that this second 
scenario is insensitive to the value of the mass m. 

We will mimic these two gedanken experiments in our analysis of open string dynamics in 
the AdSs black hole background. In Section 2 we introduce notation, describe the geometry 
explicitly, and derive the relevant equations of motion for an open string. We examine single 
quark solutions in Section 3. We first find and discuss a stationary analytic solution which 
may be viewed as describing a quark, of any mass, moving in the presence of a constant 
external electric field whose energy (and momentum) input precisely balances the energy and 

4 See Ref. [39] for an interesting qualitative attempt to understand jet quenching via AdS/CFT. 
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momentum loss due to plasma damping. Hence, this stationary solution provides the answer 
to the first gedanken experiment, and yields a measure of the (effective thermal) quark mass 
times the damping rate ji. We then turn to the analogue of the second gedanken experiment, 
and analyze the late time behavior of a moving quark in the absence of any external force. 
Looking at the low velocity, late time behavior allows us to linearize the string equation of 
motion about the static solution, reducing the problem to a calculation of the quasi-normal 
modes of the resulting linear operator. This second gedanken experiment yields the damping 
rate (i directly. 

Both the stationary analytic solution, and the quasi-normal mode analysis involve strings 
which are sensitive to the geometry arbitrarily close to the black hole horizon. This near 
horizon dependence turns up a number of subtle issues involving infrared sensitivity and 
the extent to which the total energy of a moving quark is well-defined. These issues are 
discussed in Section 3, but to insure that our interpretation is sensible, in Section 4 we 
study quark-antiquark solutions, or string configurations in which both endpoints lie on the 
D7-brane. These mesonic configurations avoid the infrared subtleties of the single quark 
solutions, but at the cost of requiring the numerical solution of nonlinear partial differential 
equations. Nevertheless, we are able to find back-to-back quark/antiquark solutions with 
external forcing in which the quark and antiquark move apart at constant velocity, as well as 
unforced back-to-back solutions in which the quark and antiquark decelerate while separating 
monotonically. The damping rate may be extracted from these quark/antiquark solutions 
when the particles are widely separated, and the results confirm and extend the previous 
analysis based on single quark solutions. 

We discuss some conceptual issues, including the effects of fluctuations which must in- 
evitably accompany dissipation in Section 5. A number of other interesting analytic solutions 
are briefly described in Appendix A. Appendix B presents the result of the quasi-normal 
mode analysis in three dimensions, where a completely analytic treatment is possible. The 
final Appendix C briefly discusses integration technique and numerical error in the non-linear 
PDE solutions of Section 4. 



1.4 Summary of results 

For a quark moving with arbitrary velocity v, we find that the rate at which it loses energy 
and momentum to the plasma is given by 



dp _ 1 dE_ 
~dl ~ v ~dt " 2 



This momentum loss rate (or viscous drag) , as a function of velocity, is independent of the 
quark mass. Note that the viscous drag may also be interpreted as the energy loss per unit 
distance traveled, since j§ = ^ 

Re-expressing the viscous drag in terms of momentum, instead of velocity, requires knowl- 
edge of the dispersion relation relating the energy E and momentum p, and hence the relation 
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between velocity v = dE/dp and momentum. As discussed below, the scale of thermal cor- 
rections to the energy of the quark in the strongly-coupled N = 4 medium is given by 



Am(T) = \\f\T . (1.3) 

A heavy quark should be viewed as one whose mass m is large compared to Am(T) (not just 
large compared to T). 5 Ifm> Aro(T), then thermal corrections to the zero-temperature rel- 
ativistic dispersion relation are negligible. In this regime, the viscous drag (1.2) is equivalent 
to 

dp 



with a friction coefficient 



Am(T) 

H = ttT — -^— L [heavy quarks, m > Am(T)\ . (1.5) 

The momentum independence of the friction coefficient is a surprising result which differs 
from the behavior of a weakly coupled plasma. 6 

The dispersion relation of lighter quarks, for which the ratio m/Am(T) is of order one, 
is substantially influenced by the medium. A quark at rest in the medium corresponds to a 
straight string stretching from the D7-brane down to the horizon. Such a quark, immersed in 
the thermal N = 4 medium at temperature T, has a rest energy M res t (T) which differs from its 
Lagrangian mass m. Determining this relation requires solving (numerically) for the change 
in the embedding of the D7-brane induced by the black hole horizon [41,42]. Asymptotically, 7 

„„, (r) = m h-^21 + U^Ilf _ JL « ) 8 + o \(^Hf] } . (, 6 ) 

[ m 8\ m / 128 V m J \_\ m J \ } 

Truncating after three terms gives a result which is accurate for M rcst + Am to within 1%. 

As m/Am(T) approaches a critical value of approximately 0.92, the thermal rest mass 
nearly vanishes. 8 Our semiclassical string analysis is only valid when the zero temperature 
mass exceeds this critical value. The resulting dependence is plotted in Figure 1. 



5 In theories with unbroken M — 2 supersymmetry the Lagrangian (or bare) mass does not get renormalized. 
So there is no need to distinguish between bare and renormalized mass, or deal with scale dependence in the 
mass — it is a physical parameter. The heavy mass regime, m Am(T), may equivalently be viewed as the 
low temperature regime, T <C 2m/\/A. In this form, one sees that the relevant scale which distinguishes low 
versus high temperature is not m, but rather m/VX This is the scale of the mass of deeply bound qq states 
which form in zero temperature M = 4 SYM with massive fundamental representation quarks [40]. So low 
temperature corresponds to the regime where these mesonic bound states form a dilute, non-relativistic gas. 

6 For a heavy quark moving through a weakly coupled plasma, the dominant mechanism of energy loss is 
two body scattering off thermal quarks and gluons, provided 1 — v 2 is not parametrically small. The resulting 
loss rate is a non-trivial function of velocity [11]. Only for small velocity is the energy loss rate well modeled 
by Eq. (1.4) with a constant value of fi. For ultra-relativistic quarks the dominant scattering process switches 
from two-to-two scattering events, in which the momentum transfer is a small fraction of the heavy quark 
momentum, to gluon bremsstrahlung in which each gluon emission can change the heavy quark momentum by 
an 0(1) amount. In this regime, characterizing the energy loss as a smooth differential process, as in Eq. (1.4), 
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1 2 3 4 5 
m/Am(T) 



Figure 1: The thermal rest mass (or energy) M rest and the kinetic mass M^ n of a quark immersed 
in the TV = 4 plasma at temperature T, as functions of the zero-temperature mass m, with all masses 
measured in units of Am(T) = ^V\T '. At m w 0.92 Ato(T), the location of the D7-brane jumps 
discontinuously to the horizon [41-44]. 



The dispersion relation of a moving quark in the thermal medium need not be Lorentz 
invariant since the plasma defines a preferred rest frame. For non-relativistic motion, the 
dispersion relation will have the form 

E(p) = M rest (T) + 2M ^ (T) + 0(p 4 ) . (1.7) 

The effective kinetic mass M^{ n (T) is not the same as the thermal rest mass M rest (T). For 
heavy quarks, we find that M k j n (T) differs negligibly from the thermal rest mass, 

/Am(T)f 



M kin (T) = M rcst (T) + o\m(-^A 

L \ m / 



but the difference between these masses becomes significant as the quark mass decreases. The 
dependence of the kinetic mass M k i n (T) on the quark mass is also shown in Figure 1. As 
ra/Am(T) approaches the lower critical value of 0.92, the kinetic mass M kui (T) has a limiting 
value just slightly greater than ^Am(T). As m/ Am(T) — > oo, both M k j n and M rest approach 
m — Am(T). 

For not-so-heavy quarks moving relativistically, we can only infer the dispersion relation 
from analysis of the time-dependent numerical solutions discussed in Section 4. We do not 
have any analytic derivation, but all our numerical results are consistent with the thermal 
dispersion relation 



E{p) = M rcst (T) - M kin (T) + yV + M kin (T) 2 , (1.9) 



no longer makes sense. 

7 The coefficient of the (Am/m) 8 term in Eq. (1.6) was determined by a numerical fit, not analytically, and 
may not be exact. 

8 The minimal value of M rcst at this point is 0.02Ara(T). When one decreases the mass beyond this point 
the location of the D7-brane jumps discontinuously to the horizon and M rcst vanishes. See Refs. [41-44] for 
discussion of this transition. 
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which reduces to Eq. (1.7) for low momentum, and gives v = dE/dp = p/\/p 2 + Mki n (T) 2 . 
For this relation between velocity and momentum, the viscous drag (1.2) is equivalent to 
p = —fip with a friction coefficient 

_ Am(T) vr ^/\T 2 
» = nT M^) = 2M^T)- (L10) 

As the quark mass decreases, the kinetic mass M^ n (T) has a lower limit of ^Am(T), and 
hence the friction coefficient has a remarkably simple upper limit, 9 

j u<2 7 rT, (1.11) 

which turns out to be dimension independent. It is tempting to speculate, along the lines of 
Ref. [34], that the ratio ji/T is bounded above by 2tt even in more general theories. 

Knowledge of the viscous drag on a quark is equivalent to knowledge of the diffusion 
constant for quark "flavor". The relation is D = T/(nM\\ n ), so our result (1.10) is equivalent 
to a flavor diffusion constant 

D = - 1 — = — jt- . (1.12) 

As discussed in Section 5, this same information may also be recast as the rate of change of 
the mean square transverse momentum of a quark initially moving in a given direction, 10 

j t {(PL?) = = 4vrAm(T)T 2 = 2vrv^T 3 . (1.13) 

This quantity, divided by the velocity of the quark (to give a rate of change per unit distance 
traveled) is sometimes called the "jet quenching parameter" q [45]. 

Physically characterizing the mechanism responsible for the energy loss (1.10) in terms of 
some microscopic picture of the dynamics of the N = 4 SYM field theory is a challenge. In the 
AdS dual description, energy and momentum flow along the string which hangs down from 
the quark, away from the D7-brane and toward the black hole horizon. It is clear that the 
portion of the string which lies close to the horizon should be thought of as describing long 
distance deformations of the medium surrounding the quark. The energy loss should not be 



9 Strictly speaking, /x equals 2nT only in the limit where M rest = 0, i.e., for an unstable D7-brane configu- 
ration. As the D7-brane is stable already for M ros t = 0.02Am(T), this limiting value of fj, is very close to the 
actual value for the lightest accessible quark masses. 

10 This result for the rate of change of mean square transverse momentum follows from modeling the effects 
of fluctuations in the momentum of the quark by a simple Langevin equation in which the noise, characterizing 
stochastic fluctuations in the force exerted on the quark, has an isotropic velocity-independent variance. At 
weak coupling [11], the friction coefficient and the stochastic force variance both show significant velocity 
dependence unless d < 1. Our strong coupling result (1.10) for the friction coefficient is velocity independent 
and valid for arbitrary values of the quark's rapidity, but we do not have a direct determination of the stochastic 
force variance for arbitrary velocities. If the force variance shows the same velocity independence as the friction 
coefficient, then the result (1.13) will also be valid for arbitrary rapidity, but if this is not true, then the result 
(1.13) will only be valid for non-relativistic motion with ti<l. 
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regarded as resulting from scattering off excitations in the thermal medium. Any scattering 
would correspond to small fluctuations in the string worldsheet, and these fluctuations are 
suppressed by inverse powers of the 't Hooft coupling (or string tension). Nor is the energy 
loss due to radiation of glueballs, which correspond to closed strings breaking off of the 
open string and whose emission is suppressed by powers of 1/N C . For relativistic velocities, 
v — ► 1, there is nothing in the classical string dynamics which is reminiscent of near-collinear 
gluon bremsstrahlung which (at weak coupling) can cause a fast moving quark to lose an 
0(1) fraction of its momentum in a single scattering. Rather, the energy transfer from the 
moving quark to the surrounding plasma should be regarded as analogous to the formation 
of a wake in the coherent polarization cloud surrounding a charged particle moving through 
a polarizable medium, or the wake on the surface of water behind a moving boat. 

Ultimately, the energy transfered to the medium from the moving quark must appear as 
heating and outward hydrodynamic flow in the non-Abelian plasma surrounding the quark. 
To see this flow directly, one would like to evaluate the expectation value of T^ v {x) in the 
presence of the moving quark. The AdS/CFT correspondence provides a recipe for this 
calculation, but its implementation is difficult. The expectation value of T^ v corresponds 
to boundary fluctuations in the gravitational metric. Deriving these fluctuations from the 
energy distribution of the string requires graviton propagators in the black hole geometry, for 
which no closed form analytic expression is currently available. Performing the computations 
required to evaluate {T^ v (x)) would allow one to examine the form of the wake behind a 
moving quark and, for example, see the sonic boom produced by a quark moving faster than 
the speed of sound. Sadly, we leave such a study for future work. 

2. Open string dynamics in the black brane background 

The AdS/CFT correspondence [30-32] posits an equivalence between A^ = 4 SU(N C ) super- 
symmetric Yang-Mills theory and type IIB string theory in a AdS^ x S 5 background. Type 
IIB strings are characterized by two numbers: a string coupling g s and a tension To, or equiv- 
alently a fundamental string length scale i s = (2ttTq)^ 1 ^ 2 . The background is characterized 
by the radius of curvature of the AdS§ and of the S 5 , which must be equal and will be denoted 
by L. 

The AdS/CFT correspondence provides a dictionary between these two seemingly very 
different physical theories. One important entry in this dictionary is the relationship between 
the string coupling and the Yang-Mills coupling, 



and an equally important entry is the relationship between the string tension, the radius of 
curvature L, and the 't Hooft coupling, 



4irg s = 9ym , 



(2.1) 



(L/4) 4 = \ = g 2 YM N c . 



(2.2) 
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AdS 


AA=4 SYM 


quantity 


L 


— 


AdSs and 5 5 curvature radius 


is 




fundamental string scale [ = Vc/] 




A 


't Hooft coupling [ = gy M N c ] 


T 




string tension [ = (2tt£ 2 )~ 1 ] 


9s 


1 2 
ii #YM 


string coupling 


Uh 


TFT 


black hole horizon radius ( x L~ 2 ) 




T 


temperature 


u m 


^ (M rcst +Am) 


minimal radius of D7-brane ( x L~ 2 ) 


T L 2 u h 


Am(T) 


thermal rest mass shift [ = ^\f\T] 


T L 2 (u m -u h ) 


M rest (T) 


static thermal quark mass 



Table 1: AdS/CFT translation table. The static thermal quark mass M rest (T) is the free energy 
of quark at rest in the TV = 4 SYM plasma. It equals the Lagrangian quark mass m in the zero 
temperature limit. 



For the convenience of readers who are not completely conversant with the AdS / CFT corre- 
spondence, these relations, plus additional ones which will appear as we progress, are sum- 
marized in Table 1. 

2.1 Adding a black hole 

The gravity dual of finite temperature J\f = 4 SU(N C ) super Yang-Mills theory is S 5 times the 
five dimensional AdS-black hole solution [33]. This solution is a geometry in which a black 
hole (or more properly, a black brane with a flat four dimensional horizon) is placed inside 
AdS space. The metric of the resulting AdS black brane solution in d + 1 dimensions may be 
written as 

ds 2 = L 2 - h{u) dt 2 + u 2 bij dxUxA , (2.3) 

\n(u) J 



where 



h(u) = u 



(2.4) 



Since the case of arbitrary dimension is usually as easy to compute as the specific d = 4 case 
of interest, we will leave d arbitrary in much of this section. Our radial coordinate u has been 
rescaled by a factor of L~ 2 ; some authors use r = L 2 u instead [46]. The black hole horizon 
is located at u = Uh where h(u) vanishes. 

The Hawking temperature of the black hole equals the temperature of the field theory 
dual. The horizon radius is related to the Hawking temperature by 



T=±« h , (2.5) 



or Uh = 7r T in d = 4. 
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Introducing a flavor of fundamental representation quarks corresponds, in the gravity 
dual of the four dimensional field theory, to the addition of a D7-brane [38]. This D7-brane 
wraps an S 3 inside the transverse S 5 and fills all of the asymptotically AdS space down to 
some minimum radial value u = u m . We require that u m > Uh- Given an open string that 
ends on this D7-brane, the quark is reinterpreted as the string's endpoint. The choice of u m 
is equivalent to the choice of mass for the quark; the relation between u m and quark mass 
will be discussed in Section 3. 



2.2 String equations of motion 

The dynamics of an open string ending on the D7-brane depends on the background geometry, 
but the back reaction of the string on the geometry is negligible and may be neglected. The 
negligible back reaction reflects the fact that fundamental representation quarks only make 
an 0(N C ) contribution to the free energy which is small, in the N c — > oo limit, relative to the 
0(N 2 ) contributions of the adjoint representation fields of JV = 4 SYM. 

The dynamics of a classical string is governed by the Nambu-Goto action, 



-T 



J da dr\f- 



detff ab . 



(2.6) 



The coordinates (a, r) parametrize the induced metric g a b on the string world-sheet. Let 
X^((j,t) be a map from the string world-sheet into space-time, and define X = d T X, X' = 
d a X, and V • W = V fM W / G fll/ where G^ v is the space-time metric (2.3). Then, writing 
det g a b = g, one has 

-g= {X-X'f- {X'f{Xf . (2.7) 

We will limit our attention to strings which lie within a three dimensional slice of the asymp- 
totically AdS space in which all but one (call it x) of the transverse coordinates x l vanish. So 
X(a, t) will be a map to (t, u, x). Choosing a static gauge where a = u and r = t, the string 
worldsheet is described by a single function x(u,t). With this choice, one finds that 



X ■ X' = L 2 (u 2 ±x f ) , 
(X'f = L 2 [h- l + u 2 {, 
{X) 2 = L 2 [-h + u 2 (x 



/\21 



(2.8a) 
(2.8b) 
(2.8c) 



and the induced metric becomes 



9ab = L 



-h + u 2 (x) 2 u 2 xx' 



u 



2 xx' \+u 2 {x') 2 



(2.9) 



The determinant of g a b is 



-JL = l-h- l u 2 (x) 2 + hu 2 (x') 2 



(2.10) 
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From this determinant, the equation of motion that follows from the Nambu-Goto action is 

x' \ u 2 d 



d_ 

du 



hu z 



-g 



h dt 



X 



-fj 



. 



(2.11) 



Useful to us in the following are general expressions for the canonical momentum densities 
associated to the string, 



7T 



TT 



-To G 



{x-x'){x v y- {x'f{x 



[IV 



-Tq G 



{X-X'){X v )-{X) 2 {X h 



[IV ' 



'-9 



(2.12a) 
(2.12b) 



For our string, these expressions reduce to 



xv? h 1 
—x x' u 2 h^ 1 
-l-(x') 2 u 2 h 



/tt! 



TqL 4 



—x' u 2 h 
-l + (x) 2 u 2 h- 1 
xx 1 u 2 h 



(2.13) 



The density of energy and x-component of momentum on the string worldsheet are given 
by 7if and 7r°, respectively. Integrating them along the string gives the total energy and 
momentum of the string, 



E 



3. Single quark solutions 



- J dorf , p = J 



da it„ 



(2.14) 



3.1 Static strings 

Single quark solutions correspond to strings which hang from the D7-brane down to the 
black hole horizon. The simplest solution to the string equation of motion (2.11) is just a 
constant, x(u,t) = xq- This solution describes a static string stretching from u = u m straight 
down to the black hole horizon at u^, and clearly represents a static quark at rest in the 
thermal medium. We may compute the energy and momentum of such a configuration using 
Eq. (2.14). The energy 

(3.1) 



E = ToL 2 / du = T L 2 (u m - u h ) , 

J Uh 



while the total momentum p (and momentum density 7r°) vanish. This energy must equal the 
(Lagrangian) mass m of the quark in the zero temperature limit. Recalling, from Eq. (2.5), 
that Uh is proportional to the temperature, we see that 



ToL u m = m 



(zero temperature) . 



(3.2) 



Moving the D7-brane to a larger radius (larger u m ) increases the mass of the quark; a D7- 
brane sitting at the boundary of the (asymptotically) AdS space corresponds to quarks of 
infinite mass. 
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However, raising the temperature affects the relation between the Lagrangian mass m 
and the position u m of the D7-brane in the gravitational description. 11 The result has the 
form 

T L 2 u m , , [T L 2 u h ^ 



i + g(^), (3.3) 



m \ m 

with the correction g(x) behaving (for d = 4) as 

g{x) = \x i -^ + 0{x 12 ). (3.4) 

Retaining just the first two terms in g(x) gives a result for T L 2 u rn /m which is accurate to 
within 1%. 

The energy (3.1) of the static string should be interpreted as the free energy of a static 
quark sitting in the thermal Af = 4 SYM medium. 12 In a small abuse of language, we will 
refer to this free energy as the static thermal mass, M rest (T). Using Uh = irT [from Eq. (2.5)] 
plus T L 2 = the result is 



M rcst (T) = T L 2 (u m -u h ) 

'Am(T)- 

m 
with 

Am(T) ee TqLV = i^AT. (3.5b) 



m-Am(T)+ mg f Am(r) ) , (3.5a) 



This relation between the static thermal mass M rcst (T) and the Lagrangian mass m is plotted 
in Figure 1. 

3.2 Moving, straight strings 

A rigidly moving string profile x(u, t) = xq + vt is also a solution to the string equation of 
motion (2.11). However, such rigid motion of the string is not physical. The problem is that 
—g is not positive definite for this profile. One finds that g vanishes at a critical value u c 
given by 



(u c r = . (3.6) 



1-v' 2 

For any non-zero velocity, u c > Uh and —g is negative in the region Uh < u < u c between 
the horizon and the critical value of the radius. A negative determinant is often a signal 
of superluminal propagation. When g = the induced metric on the string world sheet is 



n The D7-brane is a dynamical object whose equations of motion are equivalent to minimizing its worldvol- 
ume. The tension of the brane is a negligible perturbation to the background black hole geometry (suppressed 
by 1/N C ), but the D-brane does respond to variations in the background geometry by changing its embed- 
ding. According to the AdS/CFT dictionary the mass associated to a given embedding is determined by the 
asymptotic form of the D7-brane configuration. For zero temperature this procedure reproduces Eq. (3.2). 

12 This static string configuration describes the expectation value of a fundamental representation Wilson 
line in thermal M — 4 SYM [47,48]. In the finite temperature field theory, this expectation value gives e~ ,3Fq , 
where F q is the free energy, not the internal energy, of a static quark. 
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degenerate, and if —g < then the action, energy, and momentum all become complex, which 
means this solution must be discarded. By choosing x = vt, we picked inconsistent initial 
conditions where parts of the string have a velocity faster than the local speed of light at 
t = 0. While time evolution of this initial configuration gives a very simple solution, it is not 
physical. 

3.3 Moving, curved strings 

To find a physical configuration which corresponds to a quark moving at constant velocity, 
we will look for stationary solutions of the form 

x(u, t) = x(u) + vt . (3-7) 

For string profiles of this form, x', x, and yj—g are time-independent. The time derivative 
term in the equation of motion (2.11) completely drops out and we are left with an ordinary 
differential equation 

d ( , 2 x ' 



, . hu z -== = (3.8) 



where 

_/\21 1/2 



v/— 3 = L 2 [l — h- 1 u 2 v 2 + h u 2 {x'f} 1/z . (3.9) 

This differential equation is straightforward to integrate. (The v = limit of this equation 
appeared in the finite temperature calculation [47,48] of the AdS Wilson line [49,50].) The 
first integral is 

(3-10) 



or 

_'\2 



C 2 v 2 L , , ^l- 2 l-v 2 -(u h /u) d 



( x >y = ^- l-( Uh /uf j-^^--^— — 7 , (3.11) 

v 1 u s I y 1 ' J 1 - C 2 v 2 u- A - {u h /u) d ' v ; 

where C is an integration constant. This constant determines the momentum current flowing 
along the string. From the expressions (2.13) for these currents, we see that 

irl = -T L 2 C v , tt\ =T L 2 Cv 2 , (3.12) 

showing that these two currents are constant along the length of the string. 
Solving for —g yields 

-g_ = 1 - v 2 - {u h /u) d 
L 4 l-C 2 v 2 u- 4 -( y u h /u) d ' 

Both numerator and denominator are positive for large u, and negative for u near u^. So the 
only way for —g to remain positive everywhere on a string that stretches from the boundary 
to the horizon is to have both numerator and denominator change sign at the same point. 13 



(3.13) 



13 There are other solutions to these equations which have a turn-around point and do not correspond to 
strings running from large u down to the horizon. We discuss these other solutions in Appendix A. 
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Figure 2: A plot of the function F(u) which determines the string profile, in units where Uh = 1. 



This condition uniquely fixes C up to a sign, 



C = ± 



2/d 



d=4 



±- 



For the specific case of d = 4, —g/L A reduces to 1 — t> 2 , and we have 

..2 

x'(u) = ±t> ■ 



(3.14) 



h(u) u 2 

Integrating x' yields solutions of the form 

x(u, t) = x± (u, t) = xq ± v F(u; v 2 ) + vt , 



(3.15) 



(3.16) 



where the upper sign corresponds to C positive; the lower to C negative, and xq is the 
position of the string at the u = oo boundary at time zero. For arbitrary d, F(u; v 2 ) is 
a hypergeometric function. For d = 4, this function reduces to the velocity independent 
expression, 



F(u) = 



1 



2u h 



7T / U \ , / U 

arctan — — arccoth — 

2 \u h J \u h 



(3.17) 



This function is plotted in Figure 2. It vanishes as u — > oo and diverges to — oo as u — > u^; 
its asymptotic behavior is 



F(u) 



3u 3+ 

1 \n( 2Uh \ | - | 0( U ~ Uh ) 
4u h \u-u h ) 8u h \ ul ) ' 



u — > oo ; 

Uh- 



(3.18) 



The rate at which energy flows down the string is given by it}. As seen in Eq. (3.12), 
this energy flux is proportional to C. If C is positive, then energy flows down the string 
toward the horizon, and the string profile resembles a tail being dragged behind the moving 
quark, as illustrated on the left in Figure 3. If C is negative, then one has the time-reversed 
situation: energy flows upward from the horizon and the tail of string leads the quark. We 
postulate (as in Ref. [51]) that the physical process we want to describe requires us to pick 
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Figure 3: Schematic drawing of the physical solution (left) in which energy flows toward the horizon, 
and the unphysical energy solution (right) in which energy flows away from the horizon. 



purely outgoing boundary conditions at the horizon. 14 We have to discard the time-reversed 
solution in which energy flows from the black hole to the moving quark. 

The resulting rates at which energy and momentum flow toward the horizon are 



and 

— IT 



I =T L 2 u 2 h - V ^-r A ^ -VAT 2 ^=, (3.19a) 

i| =T L 2 u 2 h - V —— *3 fv^T 2 — ^— , (3.19b) 



2 



respectively. 

The stationary solution x + (u,t) given by Eqs. (3.16) and (3.17) describes an open string 
which runs from the AdS boundary at u = oo and asymptotically approaches the horizon at 
u = Uh- By truncating the solution at an arbitrary radius u m > Uh, one may equally well 
regard it as describing an open string running from a D7-brane with minimal radius u m down 
to the horizon. The rates (3.19) at which energy and momentum flow down the string are 
completely independent of u m . 

Standard Neumann boundary conditions would demand that the momentum flux 7T 1 , 
vanish at the flavor brane — so our solution with a constant non-zero it], does not satisfy 
Neumann boundary conditions. The solution is physical, but there must be a force acting on 
the string endpoint and feeding energy and momentum into the string. A constant electric field 
on the flavor brane provides precisely such a force. The field alters the Neumann boundary 
condition to the force balance condition -k\ = — F tx . 15 

14 Here and elsewhere, outgoing refers to moving out of the physical region and into the black hole. 

15 D-branes in string theory naturally support gauge fields living on their worldvolume under which the 
endpoints of strings are charged. Turning on the worldvolume gauge field will once more change the embedding 
of the D7-brane and hence the relation between u m and m. Since the results we find for the motion of the 
string in the presence of the external field is independent of u m (and hence m), the precise relation between 
the two in the presence of the field is not important for our purposes. It is, however, worthwhile noting that 
the electric field on a D-brane cannot become larger than a critical value, at which point the force pulling the 
endpoints of a string apart due to the field wins against the string tension and the system becomes unstable. 
This happens when (F tx /To) 2 = \gttQxx\ = u 4 — u\. For small mass the brane will extend down to a small 
value of u, so there will be a maximum field that the brane can support (and hence a maximal velocity) that 
decreases as the quark mass decreases. 
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Although the flux of energy and momentum in this stationary solution is finite, the total 
energy and momentum of the string is infinite, due to the contribution to the integrals (2.14) 
close to the horizon. If one simply inserts a lower limit u m \ n > Uh, together with an upper 
limit equal to the radius u m of the D7-brane, then the resulting energy and momentum (in 
d = 4) are 

/•Um ^ 

E = - du 7r t °(u) = T L 2 2 [u m - u min + t; 2 A(u m i n )] , (3.20a) 



p = du vr^(ti) = T L 2 - [u m - u min + A(u min )] , (3.20b) 

' n V 1 - V 2 



where 



A(«min) = -J- 



„ , u min „ , U m {u m ~\~Uh)(u m i n Uh) 

A arctan — A arctan — in — — — 

U h U h {U m - U h ){u min + U h ) 



(3.21) 



As n m j n — > Uh, this function diverges logarithmically as — \ Uh ln(tt m i n — Uh). The interpreta- 
tion of this IR divergence will be discussed momentarily. 

Field theory interpretation 

We want to interpret this stationary string solution as describing the steady-state behavior of 
a massive quark moving through the N = 4 plasma under the influence of a constant electric 
field 16 £ = F tx . The quark velocity will asymptotically approach an equilibrium value v at 
which the rate of momentum loss to the plasma is balanced by the driving force exerted by 
the electric field. The rate at which the electric field does work on the quark, namely v ■ £, 
should also equal the rate at which the quark loses energy to the medium. The results (3.19), 
plus the force balance condition at the string endpoint, 

irl = -S, (3.22) 

are completely consistent with this interpretation, provided one regards energy and momen- 
tum flow toward the horizon as energy and momentum transfer to the thermal medium. 

If the quark behaves as an excitation with some effective mass M and momentum p = 
Mv/Vl — v 2 , then the result (3.19b) for the momentum transfer rate is equivalent to a 
momentum loss rate dp/dt = —fip for the quark, with 

liM = f \/AT 2 = vrT Am(T) . (3.23) 

Just as in the toy model discussed in the Introduction, the momentum flow of our steady 
state solution determines /iM, but not /x or M individually. Note that \xM is independent of 
the flavor brane location u m , and hence is independent of the physical quark mass. 



16 £ is a "real" electric field — a U(l) gauge field coupled to quark flavor, having nothing to do with the 
SU(N C ) gauge fields. Such an electric field acts on fundamental representation quarks, but has no effect on 
any of the JV = 4 SYM fields. 
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The energy of the string should be regarded as the total excess free energy of the system 
- the free energy minus its equilibrium value at the given temperature. In other words, the 
energy of string includes all the effects of the disturbance to the N = 4 plasma produced by 
the moving quark. The stationary moving string solution is describing a system in which 
a quark has been forcibly dragged through the plasma (which is infinite in extent) for an 
unbounded period of time. The constant rate of work done by the external electric field thus 
translates into an infinite input of energy to the plasma. This unbounded input of energy 
is the physical origin of the IR divergence in the string energy, as may be seen explicitly 
by noting that the logarithmically divergent function A(u m i n ) appearing in the energy and 
momentum (3.20) is proportional to the difference in position (in x) between the two ends of 
the cut-off string, 

(3.24) 
with 

Ax(u min ) = x + (u m , t) - x+(u min , t) . (3.25) 

Comparison with Eq. (3.19) shows that the cut-off string energy and momentum are just a 
boosted static energy plus the net input of energy and momentum required to move the quark 
a distance Ax at velocity v, 

7-! rri j 2 ( u m — u min) 1 dE 

E = T L — - H — Ai n min , (3.26a) 

VI - v 2 v dt 

2 v(u m -u min ) 1 dp , \ 

P = T ° L /, 9 h ~ ^7 Ax ( u mm) , (3.26b) 

V 1 — v z v at 

where dE/dt = -k\ and dp/dt = —-k\ are the rates at which the external electric field transfers 
energy and momentum to the quark. Note that the u m \ n — > Uh limit of TqL 2 (u m —u m \ n ) is 
just the static rest energy M res t(T). 

The result (3.26) suggests that one might interpret M res t (T) / V 1 — v 2 as the energy (and 
v times this as the momentum) of a quark moving at velocity v through the plasma. Or 
equivalently, that the appropriate thermal dispersion relation is just a relativistic dispersion 
relation but with mass M rest (T). This, we believe, is too simplistic. A quark moving through 
the thermal plasma is a quasiparticle — an elementary excitation of the system with a finite 
thermal width given by the damping rate [i. The (free) energy of a quark moving through 
the medium is only defined to within an uncertainty given by its thermal width. A natural 
operational definition is to start with a static quark, at rest in the thermal medium, turn on an 
electric field which accelerates the quark to the desired velocity in some time r, and then define 
the energy of the quark as its initial rest energy plus the work done by the electric field while 
accelerating the quark. The acceleration time r should be small compared to the damping 
time to avoid counting energy which has already been lost to the medium. But r should 
also be large compared to the inverse kinetic energy of the quark, to minimize the quantum 
uncertainty in the energy. Satisfying both conditions is only possible if the thermal width 



ki \ 2 ^^(^min) 

A(u min ) = u h 
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is small compared to the energy, which is the basic condition defining a good quasiparticle. 
Choosing r ~ (/iMv 2 )^ 1 / 2 (for non-relativistic motion) balances the two uncertainties and 
gives a limiting precision with which one can define the kinetic energy of a moving quark that 
scales as A 1 / 4 T/p. Finding the requisite string solution with such a time-dependent electric 
field has not (yet) been done. 

3.4 Quasinormal modes 

Instead of considering a quark moving under the influence of an external electric field, we 
now turn to the motion of a quark decelerating in the thermal medium, in the absence of any 
external forcing. The setup here is the analog of the second gedanken experiment for the toy 
model discussed in the Introduction. We will focus on the late-time, and hence low-velocity, 
behavior. We extract this late-time dynamics by analyzing small perturbations to the static 
string which describes a quark at rest. A key ingredient will be the purely outgoing boundary 
conditions at the horizon which capture the dissipative nature of the process and have been 
shown to reproduce appropriate thermal physics [51]. With these boundary conditions, the 
problem becomes a quasinormal mode calculation on the string worldsheet. To complement 
the linear analysis of this section, in Section 4 we will also perform a numerical analysis of 
the full, time-dependent problem. 

The linearized equation of motion for small fluctuations around the static straight string, 
x(u,t) = xo means treating x and x' as small and retaining only terms linear in derivatives of 
x. From Eq. (2.10) one sees that this corresponds to replacing —g/L A with 1, which reduces 
the full string equation of motion (2.11) to 17 

ou v ' h 

To select the physically relevant solution, we impose purely outgoing boundary conditions 
at the horizon. These boundary conditions make the resulting boundary value problem non- 
hermitian and the resulting quasinormal modes will have real exponential time dependence. 
Close to the horizon, the most general solution of the wave equation (3.27) has the form 

x(u, t) = F(t + — log(e)) +G(t-— log(e)) . (3.28) 
V Uhd J \ Uhd J 

where e = u/uh — 1 [or u = (1 + e)uh], and F{x) and G(x) are arbitrary (differentiable) 
functions. Purely outgoing means that G = in this regime. 

Specializing to e~^* time dependence and introducing, for convenience, a dimensionless 
radial coordinate y = u/u^ and dimensionless decay rate 7 = [i/u^,, the linearized wave 
equation (3.27) becomes the eigenvalue equation 

£x = 7 2 x, (3.29a) 

with 

17 In d = 3 dimensions, this differential equation coincides with the massless scalar wave equation in the 
AdSi black hole background, at zero spatial momentum. In other dimensions, including the d — 4 case of 
interest, the linearized string equation (3.27) differs from the scalar wave equation. 
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£ = (1 _,-«) I ,.(1 (3.29b) 

Imposing outgoing boundary conditions at the horizon and Neumann boundary conditions at 
the flavor brane leads to a discrete spectrum of quasinormal mode decay rates. 

Close to the horizon at y = 1, the purely outgoing solution to Eq. (3.29) is proportional 
to (y— l) _7//rf with 7 positive. This diverges as y — > 1, showing that there are non-uniformities 
between the large time and near horizon limits. In particular, the assumption that ^J—g/L 2 ~ 
1 is only valid at sufficiently late times, when t + (u} l d)~ 1 \og(y—\) is sufficiently large. To 
evaluate the deviation of ^/^g/L 2 from unity, one needs the next-to- leading term in the near 
horizon asymptotics. Wefindx(y,t) = A(y-l)~^/ d e^ IUht {l + jB(y-l) + 0[{y-l) 2 )} with 

B = *w + a^W- This s ives 

=f = l-^ (y-l)- 2 ^ e"^ x [1 + 0{y-l)] , 

which shows that, at any fixed position outside the horizon, ^J—g/L 2 rs 1 for sufficiently large 
times. 

The near horizon asymptotics are useful also for exploring the small mass limit y m — > 1 
of our quasinormal mode problem. It is possible for x(y,t) to obey Neumann boundary 
conditions at y = y m , in the limit y m — > 1, if 7 = | + 0(y m — 1). This enables the first 
two terms in the asymptotic series expansion for x(y,t) to give comparable, and canceling, 
contributions to the slope x'(y m ,t). This result, 7 — > |, gives the intercepts of the d = 4 
and d = 2 curves shown below in Figure 4. Converting from the dimensionless decay rate 7 
back to n yields the result, valid for all d, that = 2irT in the limit where the flavor brane 
approaches the horizon. 

Large mass limit 

The differential equation (3.29) does not appear to have a simple solution for arbitrary d. 
In the special case of d = 2, the differential equation can be solved in terms of associated 
Legendre functions, as discussed in Appendix B. For the d = 4 case of interest, Eq. (3.29) is 
a particular example of the Heun equation, a differential equation with four regular singular 
points. Heun functions are difficult to work with, and other values of d appear to be even 
more difficult to treat analytically. 

In the absence of a simple analytic solution to Eq. (3.29), we attempt a power series 
solution in 7, 

x{y) = x (y) +7x1(2/) + 7 2 £ 2 (y) H . (3.30) 

We will focus on the large mass regime where the flavor brane position y m = u m /uh 3> 1, and 
we will find an iterative solution where 7 = 0(l/y m ). For the moment, the correlation between 
7 and y m may be viewed as an assumption which will be verified a posteriori. Requiring 
C x = 7 2 x implies that 

£xq = 0, £a?i = 0, and Cx2 = xq. (3.31) 
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We want to satisfy Neumann boundary conditions at the flavor brane, x'(y m ) = 0, together 
with outgoing boundary conditions at the horizon. As seen above, this requires that x(y) ~ 



A(y-l) 



-lid 



as y - 
x'(y) 



1, for some constant A, or equivalently 



A 



.1 J_ + ^My-i) +(D{ 3) 

dy-l^d? y-1 ^ 17 J 



(3.32) 



The constant xo(y) = A is the only homogeneous solution which obeys Neumann boundary 
conditions at the flavor brane at y = y m . To generate the 0(7) term in the near horizon 
behavior (3.32), the first order correction x±(y) must equal the second homogeneous solution 
of Cx = 0. The derivative of this homogeneous solution (which will be sufficient for our 
purposes) is 



x' 1 (y) = -A 



y 



d-A 



y" 



(3.33) 



At the flavor brane, x'^ym) = —Ay^ + 0(y~ 4_d ). This violates the Neumann boundary 
conditions. However, if 7 = 0(l/y m ), then this violation can, and will, be compensated by 
the next order term. Moving on to X2(y) and solving Cx2 = xq gives 



x' 2 (y) 



,d-3 



y d -l L 



^(\A^y d ) 



,d-A 



+ c 



As y 



and so 



1, the leading term in the hypergeometric function is 

^(±,1,^;/) =-i ln(i/-l)+0(l), 



= 



A 1 



In(y-l) + 0(1) 



(3.34) 



(3.35) 



(3.36) 



d 2 y-1 

The logarithmic term is precisely what is required to generate the 0(7 2 ) term in the outgoing 
boundary condition (3.32). The unwanted nondogarithmic (y— 1) _1 term is eliminated by an 
appropriate choice for the coefficient C of the homogeneous term in (3.34). 

When evaluated at y m S> 1, the hypergeometric function is negligible compared to 1, 



2 F 1 (\A^-y d m )=0{l/y m ) 



so that 



x' 2 {y m ) = Ay m 3 + 0(y m 4 ). 
The Neumann boundary condition requires 

= x'(y m ) =n/x' 1 (y m ) + 7 2 x' 2 (y m ) + 



(3.37) 
(3.38) 

(3.39) 



Inserting the explicit forms (3.33) and (3.38), one sees that the leading term from 

x'i{y m ) will cancel the leading 7 2 /y^ term from x' 2 (y m ) provided 



l=— + °Wm)- 

ym 



(3.40) 
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Figure 4: Lowest quasinormal mode decay rate 7 as a function of mass parameter y m , for d = 2 
(green dashed line) and d = 4 (thick red solid line), together with the leading large mass analytic form 
7 = — (thin blue line) 



This value of 7 2 is the smallest eigenvalue of C (for the given boundary conditions). All 
other eigenvalues are 0(1) as y m — ► 00. Expressing the result (3.40) in terms of the original 
variables, the lowest quasinormal mode decay rate is 

^ = < +0 (ul/u 2 m ). (3.41) 

The motion of the string endpoint in this quasinormal mode is a simple exponential, 

x{u m ,t) - x oc Ae~^ , (3.42) 

so the velocity of the string endpoint satisfies 

v = -pv. (3.43) 

Inserting the result (3.41) for \x and using the large mass limit of the relation (3.3) between 
the brane position and quark mass, namely u m = m/(T L 2 ), converts this result to 

^ = -(T L*ul)v, (3.44) 
with p = mv. This agrees perfectly with our earlier result (3.19b) at small v. 
Arbitrary mass 

To find the lowest eigenvalue of the quasinormal mode operator C for an arbitrary D7-brane 
location y m , it is necessary to resort to numerical analysis. A simple shooting algorithm 
suffices. At a point close to the horizon, y = 1 + e with e <C 1, one sets x(l+e) = 1 
and x'(l+e) = — This enforces the outgoing boundary condition. Then, for various 
values of 7, one integrates the differential equation (£ — 7 2 )x = out to the flavor brane, 
and successively refines 7 to locate the minimal value that satisfies the Neumann boundary 
condition x'(y m ) = 0. 
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In Figure 4 we plot the resulting lowest quasinormal decay rate, as a function of the flavor 
brane location y m , in dimensions d = 2 and 4. These numerical values are consistent with 
our large mass result that 7 = y" 1 + C(y~ 2 ). In Appendix B, we show that the d = 2 case is 
analytically soluble, and the resulting frequencies satisfy the equation 

Vm = i - \ + 0(7 2 ) • (3-45) 

When we compare the actual shape x(u,t) of the QNM with the analytic stationary 
solution, we find that they agree when 7 « 1 (that is for large mass), which is when the 
external field needed to maintain the velocity is small. 

Low velocity dispersion relation 

In addition to extracting the lowest quasinormal mode decay rate, the linearized equation 
of motion (3.27) may also be used to find the dispersion relation of a quark moving at 
low velocity. If x and x' are small, so that y 7 — g/L 4 ps 1, then the momentum density 
ir® = TqL 2 u 2 h~ l x. For a quasinormal mode with exponential time dependence, x(u,t) = 
x(u) e~ M *, the momentum density may be rewritten as 

TqL 2 2 , -1 •• TqL 2 2 /y /n A a\ 

-k x = u h x = (hu x ) , (3.46) 

where the last form follows from the linearized equation of motion (3.27). This current is 
easily integrated to find the total momentum carried by the portion of the string running 
from the flavor brane down to an IR cutoff u min > Uh, 



/Mm 
dwK° x 
-min 



p / da tt® = — — hu 2 x 



(3.47) 



Because of the Neumann boundary conditions, the upper endpoint does not contribute. Hence 

T L 2 



V 



/i 



M^min) ^min x («min) • (3-48) 



The energy may be evaluated similarly. We are interested in the deviation of the energy 
from the static rest energy, and it is therefore necessary to keep all terms up to quadratic 
order in x and x' in the energy density 7if . Suitably expanding ^—g/L 4 , the energy density 
is 

1 + \hu 2 (x'f + y 2 h^u 2 x 2 + 0(x 4 )} . (3.49) 



^0 rp j 2 

ir t = -±oL 



Using the linearized equation of motion (3.27), one may express the resulting energy solely 
as endpoint contributions, 

E = - du 7if = TqL 2 (u m - u min ) - ^h(u min ) u 2 niin x(u min ) x'(u min ) , (3.50) 
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where Neumann boundary conditions have again caused the non-static boundary term at u m 
to vanish. 
Taking 

Umin to be close to the horizon Uh, and inserting the outgoing near- horizon be- 
havior 

x{u) ~ A(u - u h )-^ dUh e-^ . (3.51) 
into expressions (3.48) and (3.50), yields a simple relation between the energy and momentum, 

p 2 

E = T L 2 (u m -u min ) + — , (3.52) 

2M kin 

where the "kinetic mass" 

T L 2 u 2 h ir^/XT 2 irT a , m , 
M kin = ^ = -1- = — Am T . 3.53 

fl 2/J, fl 

This kinetic mass is independent of the IR cutoff ii m i n , while the u m j n — > uh limit of the first 
term of the energy (3.52) is just the thermal rest energy M res t = T$L 2 (u m —Uh) ■ 

In the heavy mass limit, the decay rate \i approaches u^/u m = TrTAm(T)/m [see 
Eqs. (3.41) and (3.3)]. Therefore, in this limit the ratio of the kinetic mass M k j n (or the 
thermal rest mass M rcst ) to the Lagrangian mass m approaches one. For smaller u m , corre- 
sponding to order one values of m/Am(T), there is a more complicated relationship between 
the kinetic mass and the Lagrangian mass. The kinetic mass M kin is plotted as a function of 
m in Figure 1. 

Finally, note that the e _M * time dependence of the lowest quasinormal mode plus the 
value (3.53) for the kinetic mass are equivalent to viscous drag, 

! = -MP, (3.54) 

M = ^ = *T^. (3-55) 
2 M kin M kin V ; 

which for M = M k i n is completely consistent with the value (3.23) of \iM extracted from the 
analytic stationary solution. 

4. Quark-antiquark solutions 

Thus far, we have deduced /xM (or the viscous drag as a function of velocity) from the 
stationary analytic solution, whose analysis was valid for all velocities. And we have deduced 
the value of the friction coefficient \i itself from the linearized, quasinormal mode analysis, 
valid for late times and hence small velocities. Both approaches reveal IR sensitivity of 
the total string energy and momentum, which we argue should be viewed as reflecting an 
unavoidable level of arbitrariness in defining the partitioning of the total system energy (or 
momentum) into a piece associated with the moving quark plus a remainder associated with 
the long disturbance in the medium. 



with a friction coefficient 
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A more physical approach for dealing with this IR sensitivity is to change the question. 
Instead of considering a single quark moving through the plasma, one may study qq pair 
creation — that is, the dynamics of a quark-antiquark pair where the quark and antiquark 
are initially flying apart from each other. Such a pair corresponds to a string with both 
endpoints on the D7-brane. The previous IR sensitivity due to string dynamics arbitrarily 
close to the horizon will be cut-off by the finite quark-antiquark separation, which will limit 
how far down toward the horizon the middle of the string can "sag" . 

For simplicity, we will limit our attention to the case of back-to-back motion, so the total 
momentum will vanish and the entire string worldsheet will lie within the three-dimensional 
(t, u, x) slice of the AdS-black hole geometry. Hence, we need to find non-stationary solutions 
of the partial differential equation (2.11) with physically relevant initial conditions. To do so, 
we need time dependent numerics. 

To set up the numerical problem, it is convenient to remap the infinite range of the radial 
coordinate u £ (uh, oo) onto a finite interval. To do so, we define 



We will specialize to d = 4 and choose units where Uh = 1, so the line element of the black 
D3-brane gravitational background becomes 



where f(z) = 1 — z A and x = (x l ,x 2 ,x 3 ). Temperature dependence may be restored later 
by rescaling the coordinates: t — ► Uht, x l — > u^x 1 . In this coordinate system, the black hole 
horizon is located at z = 1 and the AdS boundary at z = 0. Our open string will end on a 
D7-brane which fills the five dimensional space from z = to z = z m . We will assume that 
both string endpoints are located at z = z m , and that the string extends only in the z and 
x = x 1 directions. 

Changing variables from u to z in the 1+1 dimensional partial differential equation (2.11), 
and then discretizing on a rectangular grid 18 in z and t turns out to be a bad approach. 
Numerical stability rapidly degrades as the string endpoints separate and the middle of the 
string gets closer to the horizon. The net result is that the numerical integration breaks down 
after a very limited amount of time. 

We have found that a much better starting point for numerical integration is the Polyakov 
action with a worldsheet metric that is a generalization of conformal gauge. Recall that the 
Polyakov action for the string takes the form 




(4.2) 




(4.3) 



As is done internally in canned PDE solvers such as Mathematical NDSolve. 



Here X(a, r) is a map from the string world-sheet into space-time, rj a p is the world-sheet 
metric, G^ v is the space-time metric, T is the string tension, and \f ZI r\ is minus the square 
root of the determinant of r) a p. We take < a < ir. 

From the action Sp, one derives the usual equations of motion for the string, 



^ d a X» d p X» ^ = d T 



G^^^X" + r ] ™X' v ) 
+ d a [G lxv ^ ] (r } ™X» + r ] ™X' v ) 

along with a constraint on the world-sheet metric, 



(4.4) 



(4.5) 



produced by the variation of Sp with respect to r] a p. The world-sheet metric rj a p may be 
integrated out (classically), converting the Polyakov action into the original Nambu-Goto 
action. 

A standard choice of world-sheet metric is "conformal gauge" , in which one chooses the 
metric to differ from a flat metric just by an overall conformal factor which is a function of a 
and r, 



\Val3\ 




j(<t,t) 



(4.6) 



Through trial and error, we have found that this choice also introduces problems for the 
numerical integration. The portion of the world-sheet close to the horizon evolves to late 
times far faster than the portion of the world-sheet closer to the boundary, introducing large 
gradients for the embedding X(a,r). A simple generalization of conformal gauge eliminates 
this problem and introduces an extra degree of freedom which may be tweaked to optimize 
the performance of the numerical integrator. Specifically, we choose a world-sheet metric of 
the form 

f -s(a,r) 



\VaP\ 







=M<t,t) 



S (7, T 



(4.7) 



We will refer to s as the stretching factor. With this choice of metric and an arbitrary 
stretching factor, the equations of motion become 



a 



A 

sz 2 
x 

sz 1 
z 

sfz 2 



da 



sft' 



( sx' 
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= o , 



sz' 
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"2s" 



(i 2 - s 2 t' 2 ) d. 
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)d z 



- (x 2 -s 2 x' 2 )d z (^ 
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(4.8a) 
(4.8b) 

(4.8c) 
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and the constraints, written explicitly, are 



= -/ it' + xx' + f" 1 zz' , (4.9a) 
= _/ (i 2 + s 2 t' 2 ) + (x 2 + s 2 x' 2 ) + r\z 2 + s 2 z' 2 ) . (4.9b) 

Here (and below), x = d T x and x' = d a x, etc. 

The next step in setting up the numerical integration is to find good initial conditions, a 
task which is made harder by the constraint equations. We have found two consistent sets of 
useful initial conditions, both inspired by the classical limit of the leading Regge trajectory 
of the string quantized in flat space. Recall that the leading Regge trajectory for an open 
string in flat space with pure Neumann boundary conditions has a classical limit, 

t = At , x = Acosa sinr , z = Acosa cost , (4-10) 

which corresponds to a line segment of length A spinning in a circle in the xz-plane. Once we 
introduce a D7-brane along z = zo, there is a closely related classical string state with mixed 
Neumann-Dirichlet boundary conditions which serves as our inspiration for initial conditions, 

t = At , x = ^cosct sinr , z = zq + A sin a sinr . (4-11) 

This describes a semi-circle expanding and contracting in the xz-plane. 

Our first set of initial conditions for the AdS-black brane geometry can be thought of as 
the r = limit of the semi-circle solution with some scaling factors that compensate for the 
fact that we are no longer in flat space, 

i(a,0) = 0, i(a,0) = Ail-z^}- 1 / 2 , (4.12a) 
x(a,0)=0, x(a, 0) = A cos a , (4.12b) 
z(a,0)=z m , z(a, 0) = A [l-z^} 1/2 sin a. (4.12c) 

We will call these boundary conditions "point-like"; at time t = 0, the string is mapped 
onto a single point in space-time. The parameter z m , which controls the AdS radius of the 
string endpoints, is determined by the quark mass, z m = Uh/u m . Notice that the initial 
speed of the ends of the string (given by x/i) is constrained to be yjl — z^, so the initial 
speed is an increasing function of quark mass. The "amplitude" A controls how much energy is 
contained in the initial (zero length) string. [The explicit relation is given below in Eq. (4.16).] 
Physically, these initial conditions should resemble the effect of a local current which produces 
a quark-antiquark pair when acting on the thermal equilibrium state, with the quarks having 
sufficient energy so that their dynamics may be regarded as classical. 
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Our second set of initial conditions is an expanding semi-circle characterized by an ad- 
justable speed v. We take 



t(<r,0) = 0, (4.13a) 

x{a,0) = Acosa, (4.13b) 

z(a,0) = z m + A sin a, (4.13c) 

x(a, 0) = v cos a i(a, 0) . (4.13d) 

The constraint equations then force 

z(a, 0) = v f sin a i(a,0) , (4.13e) 

and 

A s [sin 2 a + / _1 cos 2 a] 1 " 
[/ - v 2 (cos 2 a + / sin 2 a)] 



°) = - H o lV 2 ' (413f) 



with / evaluated at z(a, 0). For these initial conditions to result in a real valued X(a, r), the 
inequality 

v 2 (cos 2 a + fsin 2 a) < f (4.14) 

must be satisfied for all a. This macroscopic "semi-circle" initial configuration does not 
correspond to the action of any local operator, but the finite size of the string allows more 
freedom in choosing the initial speed v of the quarks. 

We used the NDSolve routine in Mathematica for numerical integration [52]. See Ap- 
pendix C for a discussion of numerical error. 

4.1 Forced motion 

We would like to confirm that the analytic solution presented in Section 3.3 is physically 
relevant. One might worry because of the IR divergence in the string energy. Therefore, we 
will investigate a quark-antiquark pair in the presence of a constant electric field, which will 
naturally drive the quark and antiquark in opposite directions. 

A straightforward numerical integration of the equations of motion (4.8) cannot handle 
sending the flavor brane all the way to z m = (or the quark mass to infinity) due to the 
divergence of l/f at the boundary. So we will choose positive values of z m and use either 
the point-like initial conditions (4.12) or the semi-circle initial conditions (4.13) to create 
a separating quark-antiquark pair at time t = 0. Instead of the usual Neumann boundary 
conditions at z = z m , we simply fix the endpoint velocity, x/t = v at a = and —v and 
a = it. This corresponds to a time-dependent electric field which is asymptotically constant, 
and whose strength is adjusted precisely to cancel the viscous drag on the quarks at all times. 
As noted above, the speed v must equal — zf n for point-like initial conditions, or satisfy 
condition (4.14) for semi-circle initial conditions. In either case, we find numerical solutions 
which nicely match onto two copies of the analytic solution (3.17) at late times. As the 
ends of the string separate, the middle of the string droops toward the horizon and the time 
dependence approaches the stationary form (3.7). 



- 28 - 




Figure 5: Cross sections of the string world sheet at times (a) t — 0.25, (b) t — 0.5, (c) t = 1.1 and 
(d) t = 1.8 (in units where Uh = ttT = 1), for point- like initial conditions, constant velocity boundary 
conditions, and a D7-brane at z m — 0.75. The velocity v/c = 0.83 and the value of z m corresponds to 
m/Am(T) — 1.28. The red dashed line shows the constant velocity analytic single quark solution. 



In Figure 5, we plot numeric results for point-like initial conditions with a D7-brane at 
z m = 0.75 and A = 0.25. This value of z m corresponds to a mass ratio m/ Am(T) = 1.28 and 
a speed v/c = 0.83. As time goes by, the expanding string, plotted in black, matches onto the 
analytic solution, plotted in red, more and more closely. By t = 1.8, the two are practically 
indistinguishable for x > 0. 

Through experimentation, we found that a stretching function of 

' (4.15) 



1 - 3.5 z(l-z) 

works particularly well for this string. The 1—z factor in the numerator prevents the string 
worldsheet from being dragged to late times close to the horizon. The denominator amelio- 
rates some distortion of the constant a and r contours of the worldsheet for intermediate 
values of z. 



4.2 Unforced motion 

In this subsection we consider the motion of a quark-antiquark pair created with point-like 
initial conditions (4.12) in the absence of any external forcing. That is, we impose the usual 
Neumann boundary conditions along the D7-brane. A short calculation shows that the energy 
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Figure 6: String endpoint as a function of time for an oscillatory solution with z m = 0.50, E/M rcst = 
1.2, and pointlike initial conditions. To observe many oscillations, this integration tolerated a larger 
numerical error (10 -3 instead of 10 -4 typical of the other plots in this section). The time dependent 
damping is comparable to the numerical error. 

[given by Eq. (2.14)] of our point-like initial configuration is 19 



The subsequent motion is sensitive to the value of this initial energy of the qq pair. The 
static potential between a heavy quark and antiquark in M = 4 SYM rises linearly at short 
distance before switching to a Coulombic form at a cross-over distance set by the inverse quark 
mass [40]. At non-zero temperature, the potential remains approximately Coulombic until a 
distance of order the inverse temperature, where (at -/V c = oo) there is an abrupt transition to 
to a constant limiting value [which is twice the static thermal rest energy M res t(T)] [47, 48]. 20 If 
the energy of the qq pair is sufficiently low, then the attractive force between the separating 
quarks will be strong enough to cause their trajectories to turn around, and the resulting 
motion will resemble an oscillator. If the energy is sufficiently high, then the quark trajectories 
will not have turning points and the motion will resemble a highly overdamped oscillator. 

For our point-like initial conditions, the exact energy threshold for non-oscillatory motion 
need not equal 2M rest (T) precisely, because of the possibility of exciting internal string degrees 
of freedom. Numerically, however, there does indeed appear to be a divergence in the period 
at E > 2M rest (T) for these point-like initial conditions. 

Numerically we do find oscillating solutions similar to the flat space expanding and con- 
tracting semi-circle (4.11) when E is well below 2M rest (T). An example is illustrated in 
Figure 6. As E approaches 2M res t(T) from below, the period of oscillation becomes longer 
and longer. 

19 The current 7r° can be obtained from a limit of (2.12a), or more directly from the Polyakov action. One 
finds tt? = -T L 2 V TT v^jG tt i . 



This corresponds to a change in the lowest energy string configuration from one in which a string connects 
the two quarks, to one in which a string runs straight down from each quark to the horizon. 



E = T L 




(4.16) 
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Figure 7: Plots of the position of the endpoint of the string versus time for several different D7-brane 
positions z m . The string endpoint is plotted in black while the best-fit curve, as described in the text, 
is a green dashed line. The two curves are nearly identical. Parameters for the different solutions are: 
(a) z m = 0.75, E/M reBt = 14, (b) z m = 0.50, E/M lcst = 6.1, (c) z m = 0.25, E/M vcst = 4.4, and (d) 
z m = 0.25, E/M Tes t = 5.9. Units where Uh = i^T — 1 are used. 



To extract information about the viscous damping of a single quark, we want to create a 
quark-antiquark pair with an energy greatly exceeding the binding energy, 

E » 2M rcst (T) , (4.17) 

so as to minimize the interaction between the quark and antiquark. Numerical solutions 
satisfying this large energy condition do show the expected non-oscillatory behavior. Several 
examples are shown in Figure 7. Solutions (a) and (b) used a stretching factor s = (1— z), 
amplitude A = 0.25, and flavor brane positions z m = 0.75 and 0.5 respectively. The masses 
and energies were ra/Am(T) = 1.28, M rest (T)/m = 0.26, and E/M rest (T) = 14 for (a) and 
m/Am(T) = 1.98, M rest (T)/m = 0.51, and E/M Iest (T) = 6.1 for (b). For (c), a stretching 
factor s = (1— z) 3 / 2 , amplitude A = 0.17, and flavor brane position z m = 0.25 were used. 
For this run, m/Am(T) = 4.0, M rcst (T)/m = 0.75, and E/M Ie8t (T) = 4.4. In the last run 
(d), a stretching factor of s = (1— z) 2 , an amplitude A = 0.20, and a flavor brane position of 
z m = 0.25 were used, corresponding to E/M rest (T) = 5.9. Numerical error limits how far we 
are able to integrate in time. For initial conditions corresponding to lighter or less energetic 
quarks, one sees a large decrease in velocity and a clear approach to an asymptotically constant 
position. But for higher energies or more massive quarks (which experience less damping) 
numerical error prevents us from following the quark to the non-relativistic regime. 
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Given a quark-antiquark pair with large energy, we model each quark independently as 
a particle experiencing a damping force 

< 418 > 

If momentum is proportional to velocity, as with non-relativistic motion, then this equation 
integrates to 

v{t) = VQe-" t , x{t) = Xoo - - e""* . (4.19) 



But if the relation between velocity and momentum has a relativistic form, p oc v/y/l — v 2 , 
then equation (4.18) is equivalent to 



^ = -»v(l-v 2 ), (4.20) 



which integrates to 



and 



v(t) =vq/ yjvl + (l-vl)eW, (4.21a) 

i + v(ty 



x(t) = Xoo - — In 



v(t) 



(4.21b) 



To test the validity of this description and extract information on the damping rate fi we fit 
the numerical results for the position of the string endpoint to Eq. (4.21), treating fi, ln(l— Vq) 
and Xoo as free parameters. 21 As shown in Figure 7, when the large energy condition (4.17) 
holds the resulting fit, using the form (4.21), is quite good for all r, although we do see some 
minor deviations at small r. In contrast, fits of these large energy solutions to the simple 
exponential form (4.19) are somewhat worse. 

The fact that fits using Eq. (4.21) are so good supports the presumed relativistic relation 
between p and v together with a friction coefficient fx that is independent of p. We are 
following the evolution of the quark over typically a rather large momentum range as can 
be seen from Figure 7. The results for the extracted values of \x are nearly independent of 
the energy so long as E > 2M rcst (T). Changing the energy by 50% or more typically only 
results in a few percent change in the value of fx. For example, for a D7-brane at z m = 0.5, 
as we increase E/M Test from 4.8 to 7.2, fi changes from 0.79 to 0.80. For the case z m = 0.25, 
changing E/M TCSt from 4.4 to 5.9 changes fx from 0.325 to 0.326. 

Table 2 compares the values of fx from the quasinormal mode calculation of Section 3.4 
to the best fits of our numerical integrations. The numbers are astonishingly close, giving us 
confidence in the linear analysis. The results begin to differ by a few percent for small mass 
quarks. It is not a-priori clear what causes this discrepancy. The extraction of the energy 



21 Given a numerical integration up to a worldsheet time T ma x, we fit the endpoint of the string to the 
assumed form for x(t) only for r mal /2 < r < r max , choosing ten equally spaced points in this region. We limit 
the data used in the fit to the latter half of the available time interval in order to minimize the effects of the 
interaction between the quarks, which is largest when the quarks are close together. 
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Table 2: The friction coefficient /j, for various values of quark mass m. The quasinormal mode (QNM) 
results were calculated using the linear analysis of Section 3.4, while the numeric results come from 
the full time dependent numerics discussed in this Section. 

loss rate may be affected by the interaction with the other quark, but that should produce 
a correction of the opposite sign. There may be small thermal corrections to the relation 
between momentum and velocity, or residual momentum dependence in the damping rate. 
Conceivably, there could be nonlinear effects that are absent in the quasinormal mode analysis 
but which reappear in this full dynamical simulation and are relevant even at (accessibly) 
late times. Or this small discrepancy may reflect residual errors in our numerical integration. 

5. Discussion 

Let us close with a discussion of the validity of our approximations. The classical treatment 
of the string is justified as long as the string is much longer than a string length; quantum 
fluctuations will be suppressed by powers of i s /R where R denotes the characteristic length 
of the string. All the single quark solutions we considered had strings with length of order 
L, the AdS curvature radius, or larger. As one lowers the quark mass toward the critical 
value m c « 0.92 Am(T), the D7-brane approaches but does not quite reach the horizon. The 
shortest string one can get (with m > m c ) has a length of about 0.02L. Since L = \i £ 8 , 
we see that for sufficiently large A, quantum fluctuations of the string are always suppressed. 
Hence, a classical treatment of the string dynamics is valid for sufficiently large A as long as 
the quark mass exceeds the critical value m c . 

For applications to QCD, however, one may be interested in large but not asymptotically 
huge values of the 't Hooft coupling, perhaps A « 20 (corresponding to a s « 0.5). In this 
regime, the condition R 3> i s can become non-trivial and for masses near the critical value 
m c (which is about 2.2 T for A « 20) quantum fluctuations of the string will be important. 

Another phenomenon that has not appeared in our discussion up to now is Brown- 
ian motion. Any dissipative thermal system must also have fluctuations, as shown by the 
fluctuation-dissipation theorem. In particular, a quark (of finite mass) initially at rest in the 
plasma should not be able to remain at rest. It will undergo Brownian motion and diffuse 



- 33 - 



away from its starting point. Over a time t it will travel a distance Ax ~ y 'Dt. The diffusion 
constant D is directly related to the viscous drag, 



This physics is missing in our classical treatment of string dynamics. A motionless string 
stretching from the D7-brane to the horizon is a solution to the equations of motion, and 
there is no obvious reason for the string endpoint to move at all. This straight motionless 
string clearly represents a quark at rest in the plasma. 

The reason we do not see Brownian motion is due to the non-uniform nature of the large 
A and large time limits. For a quark initially moving with some O(l) velocity v, stochastic 
Brownian motion will be unimportant until sufficiently late times. To see this explicitly, one 
can use the Langevin equation 

f=-»P + £(t) (5.2) 

to model the behavior of the quark. The required assumptions underlying this description 
are discussed below. Here £(t) is stochastic white noise, with variance 

m)Ut'))=CS ij S(t-t'). (5.3) 

Calculating the mean square momentum at time t gives 



(pi(t) Pj (t)) = (piWXpjit)) + ^ S i3 (1 - e- 2 ^) . (5.4) 



with (p(t)) = p(0) . At equilibrium, by equipartition, the kinetic energy p 2 /(2Mk in ) of the 
quark must be §T. (Note that the condition m>m c automatically implies that Mk in » T, so 
a non-relativistic form of kinetic energy is appropriate for quarks in equilibrium.) Requiring 
that the large t limit of (p(t) 2 } equal 3TMk; n shows that the strength of the noise is determined 
by the viscous drag, 

C = 27>M kin (T) = 2vrT 2 Am(T) = ttv^T 3 . (5.5) 

Integrating the Langevin equation (5.2) again to find the quark's position, assuming non- 
relativistic motion and vanishing initial velocity, and computing the mean square displacement 
gives 

(A Xi (t)A Xj (t)) = ^ 6ii[l + 0(l/(»t))] . (5.6) 

This must equal the classic result (Axj(t) Axj(t)) = 2Dt5ij for the variance of the proba- 
bility distribution P(Ax,t) = (4:irDt)~ 3 ^ 2 e~ Ax /( 4£) *) generated by a diffusion equation with 
diffusion constant D. Combining Eqs. (5.5) and (5.6) gives the stated value (5.1) for the 
diffusion constant. 
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For diffusive effects to be negligible, we require the second term of (5.4) to be small 
compared to the first term, giving an upper bound on the time over which a deterministic 
treatment of quark motion is valid, 

1 / 2K.E.\ , x 

4<tBS _ ln ^ 1 + _), (5.7) 

where K.E. = ^M^[ n Vq is the kinetic energy of the quark at time zero. At the same time, 
for Eq. (5.2) to adequately model the energy loss of the quark, quantum uncertainty in its 
kinetic energy must be negligible compared to the change in the kinetic energy which we wish 
to resolve. This imposes a lower limit on the time during which our classical description is 
valid, namely 1/t <C fit (K.E.), or 

t » t Q = [fi x {K.E.)]- 1 ' 2 . (5.8) 

Our analysis requires that tQ <C t <C £#. Using /j, ~ Am(T)T / M kin , this condition may 

be rewritten as 

Am /K.E. / 2K.E.\ 

«i — In 1 + — . (5.9) 



M kin V T V 3T 

This inequality is most stringent for the lightest (accessible) quarks with m near m c and 
-^kin ~ fAm(T). Hence the required large separation between the quantum and diffusive 
time scales is valid as long as the initial kinetic energy of the quark is large compared to the 
temperature, 

K.E. = ±M kin v 2 Q > T . (5.10) 

For a quark mass near the critical limit m c , this is equivalent to the requirement that the 
initial velocity satisfy 

v » X~ l/A . (5.11) 

A point to be emphasized is that the classical treatment of quark dynamics underlying all 
results in Table 2 is valid provided A is sufficiently large. 

The result (5.4) for momentum fluctuations also allows one to relate the rate of change 
of the mean square transverse momentum to the noise, and hence to the viscous drag, 

j t (P±(tf) \ t=Q = 2C = 4T/xM kin (T) = 2ttV\T 3 . (5.12) 

This characterizes the diffusion in transverse momentum of a quark propagating through the 
plasma. However, the simple connection (5.12) between transverse momentum diffusion and 
viscous drag relies on the validity of the Langevin equation (5.2) with isotropic white noise 
(5.3) for modeling the stochastic force fluctuations acting on the quark. If the stochastic force 
fluctuations have significant momentum dependence, or non-Gaussian correlations, then this 
simple description will not be adequate. One can argue that the simple Langevin description 
is valid for non-relativistic motion, »<1. Whether it remains valid for arbitrary momentum 
is not completely clear. At weak coupling, both the viscous drag and the noise variance 
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acquire significant velocity dependence for 0(1) values of rapidity [11]. Since we find no 
velocity dependence in the friction coefficient n at strong coupling, it seems plausible that 
there will also be negligible velocity dependence in the variance of the force fluctuations, even 
for relativistic motion. However, this has not been directly verified. 

Note Added: 

Several related papers [53-55] have very recently appeared which overlap with portions 
of our analysis. The quark diffusion constant found in Rcf. [54] agrees with our value (1.12). 
The result of Ref. [53] for the jet quenching parameter q does not agree with our result (1.13). 
However, these authors are addressing a different physical question involving radiative energy 
loss of a lightlike projectile. Interestingly, their result has the same parametric dependence 
as our result (5.12), but with a coefficient which is 15% smaller. 
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A. Other solutions 

In this appendix we briefly discuss other stationary string solutions in the AtiS^-black hole 
background. As we showed in Section 3.3, the ansatz x(u,t) = x(u) + vt reduces the string 
equation of motion to the ordinary differential equation (3.8) with a first integral (3.11). The 
numerator in the expression (3.11) vanishes at u c = (1 — v 2 )~ 1 / d iif l , while (for d = 4) the 
denominator vanishes at no = (u\ + C 2 v 2 ) l / A . The previously discussed single quark solution 
requires u c = no, so that x' is non- vanishing and non-singular everywhere between the horizon 
and the boundary. But if these two radii do not coincide, one can still find solutions with —g 
everywhere positive on the worldsheet. If no > n c , then the physical solution lives entirely in 
the u > «o part of space, depicted on the left in Figure 8. This configuration corresponds to 
an infinitely heavy external quark/antiquark pair at finite temperature, moving at a common 
velocity v. For a static quark/antiquark pair, it was found in Ref. [47,48] that the solution 
only exists for a bounded range of quark/antiquark separations I. At zero velocity the largest 
separation l max is achieved when no approaches Uh- At non-zero velocity, this solution ceases 
to exist beyond uq = u c . 

Alternatively, if no < n c then the solution lives entirely in the region between no and the 
horizon, as depicted on the right in Figure 8. One might think that this string solution could 
represent some coherent gluonic excitation. But since the momentum is outgoing on one end 
of the string, and ingoing on the other, we believe this solution is unphysical. 

Finally, yet another very simple solution is a straight string, x = vt, moving at constant 
velocity. We argued in Section 3 that such a string, stretching from u m down to the horizon 
and moving at any non-zero velocity, is not physical. However, if one considers a theory 
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Figure 8: Left: Semiclassical solutions corresponding to a quark/ antiquark pair with a fixed spatial 
separation moving through the finite temperature medium at constant speed. Right: A stationary 
solution in which the string moves at constant velocity outside the horizon. One end of the string 
satisfies physical (outgoing) boundary conditions at the horizon, but the other end does not. Hence, 
this solution is unphysical. 




u h 

Figure 9: A straight, moving string solution corresponding to a light-heavy meson (in a multi-flavor 
theory) moving through the thermal medium at constant velocity. 

with two flavors of quarks with different masses, so that the gravitational description has two 
D7-branes at differing radial positions u mi and u m2 , then one can regard the portion of this 
trivial solution lying between u mi and u m2 as describing a moving light-heavy meson. Figure 9 
depicts this configuration. The solution is physical provided u mi and both greater 

than u c . This condition shows that, for given quark masses, there is a maximum velocity 
with which such a color-singlet meson can move through the medium without experiencing 
any drag (at leading order in l/N c and A — > oo), 

v 2 < 1 - (u h /u m< ) d , (A.l) 

where u m . is the lesser of u mi and u mo . 



-37- 



B. Quasinormal modes in d = 2 



In d = 2, the metric function h(u) = u 2 — and the resulting linear equation (3.27) can be 
solved analytically. The most general solution is 

x(y,t)=e^ t f(y)=e-^(c P P?(y)/y + c Q Ql(y)/y), (B.l) 

where we have introduced the dimensionless quantities 7 = fi/uh and y = u/uh, and P^ 
and Qj are associated Legendre functions. (We follow the conventions in Gradshteyn and 
Ryzhik [56].) The horizon is at y = 1, and the boundary is out at y — ► 00. In order to study 
the behavior of P± and Q\ at large y, and for y close to 1, it is convenient to use the relation 
to hypergeometric functions. The associated Legendre functions P± and Qj are uniquely 
defined for 1 < y < 2 and y > 1, respectively, as 

Q?(y) = T(2+ 7 ) (l - ^ 2 ) 7/2 <T 2 2 F X (i + §, | + 1, §; y" 2 ) . (B.3) 

These expressions may be used to evaluate P^(y) for y close to 1 and Qj(y) at large y. To 
evaluate Pi(y) at large y and Qj(y) close to 1 one can use hypergeometric identities to find 
(see also Ref. [56]): 

r?(v) = ^73^ (! - ^" 2 ) 7/2 ^ 2 ^ (i+§> i+i l^" 2 ) 

7/2 

QT(y) = ^ 



r(7)(^) 7 2 F 1 (-l,2,l- 7 ;I-l) 

r(- 7 )r(2+7) /y-iy /2 r t y) 

+ r( 2- 7 ) UTiJ 2Fl t" 1 ' 2 ' 1+7 ' 3-3) 



(B.5) 



Close to the horizon we find that Q\ is a linear combination of a solution that goes as 
(y— I) 7 / 2 and one that goes as (y— I) -7 / 2 . According to Eq. (3.28), when combined with e~^ 1 
time dependence the former is in-going while the latter is outgoing. On the other hand, P± 
only has a (y— 1)~ 7 / 2 term and hence is the purely outgoing solution. So to find quasinormal 
modes we can focus on the P^ solution only, 22 and set cq = 0. 

Imposing Neumann boundary conditions at a flavor brane, the quasinormal modes are 
given by solutions to 

/^U . = W(y)/y)\ y = ym = o- (b.6) 

To check the results we also looked at e^* time dependence instead of e _Mt . In that case {y-l) l/2 is the 
physical near horizon behavior. This implies a particular linear combination of PJ 1 and QJ. The final answer 
turns out to be the same. 
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Figure 10: Quasinormal mode wavefunction f(y) for a mass corresponding to y m = 10. 



For large y m we can expand the hypergeometric functions to obtain 



r(2- 7 ) 



[i - ± 7 V 2 + o ( 7 V 4 )] 



3r(-i- 7 ) 



T -^[l + ^(3 + 9 7 + 7 2 )y- 2 + 0( y - 4 )]. 

(B.7) 

Assuming that 7 ~ we have kept all terms up to order y~ 7 . This assumption will be 

justified presently. From this asymptotic expansion, it follows that 



1X-1-7) 



I-7 2 



1 - (1 -+^)y m 3 + o(y- 5 ) 



(Bi 



Solving f'(y m ) = for y m yields the asymptotic expression y m = ^ — ^ + 0(-y 2 ), which justifies 
a posteriori our assumption about the scaling behavior of 7. Inverting this expression for 7 
yields 

7 = ^--i + Ofc 4 )- (B.9) 
The corresponding quasinormal mode wavefunction f{y) is plotted in Figure 10 for y m = 10. 



C. Numerical error 

To perform the numerical integration in Section 4, we used the NDSolve routine provided by 
Mathematica [52] and checked for numerical error in a variety of ways. We used the spatial 
error estimate provided by NDSolve, we checked the constraint equations, and we compared 
the numerical integration results for different grid spacings. 

The routine NDSolve produces a warning if its internal spatial error estimate exceeds 
a threshold. None of the numerical results we present here generated such warnings. This 
spatial error estimate is formed by considering the final step in the numerical propagation. 
One additional time step is made both with the original grid and with a coarser grid with 
half the number of grid points. Using the Richardson extrapolation formula, Mathematica 
produces a warning if the difference ^"j^ > 10, where 2/2 and y\ are the resulting values 
of the function for the two different choices of grid and p is the order of the discretization 
routine that converts derivatives into differences, the default value of which is p = 4. The norm 
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1 2/2 — 2/1 1 involves a scaled sum over the difference between yi and y\ which the documentation 
of NDSolve does not describe in detail. 

We made our own rough estimates of the numerical error by numerically integrating each 
example twice with two different grids. If the first grid was specified to have a minimum 
number of points n, the second grid would have a minimum of In points. The (absolute) 
difference between the two integrations was kept under about 10~ 4 . 

The constraint equations (4.9), if satisfied by the initial conditions, should remain satisfied 
at all later times. But unless special steps are taken when converting derivatives into finite 
differences, the differencing scheme will no longer preserve the constraint equations exactly. 
The extent to which the constraint equations are satisfied at later times is thus an indirect 
measure of the accumulated numerical error. In the numerical results presented here, we 
checked the constraints at a handful of points and found that they differed from zero only by 
about 1(T 6 . 
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